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ABSTRACT 

Context. Observations at 0/ 1 have revealed the existence of dark cores in the bright filaments of sunspot penumbrae. Expectations are 
high that such dark-cored filaments are the basic building blocks of the penumbra, but their nature remains unknown. 
Aims. We investigate the origin of dark cores in penumbral filaments and the surplus brightness of the penumbra. To that end we use 
an uncombed penumbral model. 

Methods. The 2D stationary heat transfer equation is solved in a stratified atmosphere consisting of nearly horizontal magnetic flux 
tubes embedded in a stronger and more vertical field. The tubes carry an Evershed flow of hot plasma. 

Results. This model produces bright filaments with dark cores as a consequence of the higher density of the plasma inside the tubes, 
which shifts the surface of optical depth unity toward higher (cooler) layers. Our calculations suggest that the surplus brightness of 
the penumbra is a natural consequence of the Evershed flow, and that magnetic flux tubes about 250 km in diameter can explain the 
morphology of sunspot penumbrae. 
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1. Introduction 

At high angular resolution, penumbral filaments are observed 
to consist of a central dark lane and two lateral brightenings 
(Scharmer et al. 2002; Sutterlin et al. 2004; Rouppe van der 
] Voort et al. 2004; Bellot Rubio et al. 2005; Langhans et al. 2007). 
The common occurrence of dark-cored filaments and the fact 
that their various parts show a coherent behavior have raised ex- 
pectations that they could be the fundamental constituents of the 
penumbra. Their nature, however, remains enigmatic. 

One possibility is that dark-cored filaments represent mag- 
netic flux tubes carrying a hot flow. This would support the 
uncombed model proposed by Solanki & Montavon (1993), 
which describes the penumbra as a collection of nearly hori- 
zontal flux tubes embedded in a more vertical background field. 
The uncombed model is, by far, the most successful represen- 
tation of the fine structure of the penumbra currently available. 
It explains the polarization profiles of visible and near-infrared 
lines observed in sunspots (e.g., Beck 2008), including their net 
circular polarization (NCP). This success is not trivial, since 
the behavior of the NCP depends on the details of the mag- 
netic and velocity fields in a very subtle way (see Mtiller et 
al. 2002; Borrero et al. 2007; Tritschler et al. 2007, and refer- 
ences therein). The uncombed model is supported not only by 
observations, but also by theoretical work. Schlichenmaeier et 
al. (1998) and Schlichenmaier (2002) performed numerical sim- 
ulations of penumbral flux tubes in the thin tube approximation. 
These calculations show filaments whose morphology and dy- 
namics are very similar to those actually observed in the penum- 
bra (Schlichenmaier 2003). Moreover, the simulations offer a 
natural explanation for the Evershed flow, the most conspicuous 



dynamical phenomenon of sunspots. In spite of these achieve- 
ments, it is still not known whether magnetic flux tubes can 
also account for the existence of dark cores in penumbral fila- 
ments and, more importantly, for the surplus brightness of the 
penumbra. Schlichenmaier & Solanki (2003) suggested that hot 
upflows along magnetic flux tubes would indeed be able to heat 
the penumbra to the required degree if the tubes return to the 
solar interior after they have released their energy in the photo- 
sphere. At that time the existence of opposite-polarity field lines 
in the penumbra was unclear, but now it is a well-established ob- 
servational fact: submerging flux tubes have been detected from 
Stokes inversions and even imaged directly by Hinode (Sainz 
Dalda & Bellot Rubio 2008). Thus, the Evershed flow remains 
the best candidate to explain the brightness of the penumbra. 

Another possibility is that the dark-cored filaments are the 
manifestation of field-free gaps that pierce the sunspot mag- 
netic field from below. The concept of a gappy penumbra was 
proposed by Spruit & Scharmer (2006) and Scharmer & Spruit 
(2006) as an alternative way to explain the surplus brightness of 
the penumbra, on the assumption that the Evershed flow is not 
sufficient. The gaps would sustain normal convection, carrying 
heat to the solar surface. Radiative transfer calculations need to 
be performed to show that a penumbra consisting of field-free 
gaps is able to explain the corpus of spectropolarimetric obser- 
vations accumulated over the years. In its present form, however, 
the gappy model is bound to experience substantial difficulties 
when confronted with the observations (Bellot Rubio 2007). 

Recently, Heinemann et al. (2007) have presented first at- 
tempts to simulate the penumbra in 3D. The parameters govern- 
ing the calculations are still far from those of the real sun and, 
as a consequence, the model sunspot does not show a typical 
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Fig. 1. Top left: magnetic field lines in the xz-plane. The circle centered at (0, 0) km represents the flux tube's boundary. Note the 
wrapping of the field lines around the tube. Top right: field strength distribution. Bottom left: field inclination distribution. Bottom 
right: gas density distribution, for temperatures in the tube and background given by the cool model of Collados et al. (1994). Shown 
are density differences with respect to the unperturbed atmosphere (Fig. [2]). 



penumbral pattern. Yet, an interesting result of the simulations 
is the existence of small blobs of plasma with weaker and more 
inclined fields than their surroundings. The magnetic properties 
of these structures are reminiscent of those of the flux tubes pos- 
tulated by the uncombed model. Some of them show a dark lane 
similar to the dark cores of penumbral filaments. The dark lanes 
are produced by locally enhanced density and pressure that shift 
the r = 1 level to higher photospheric layers, where the tem- 
perature is lower. This effect was identified for the first time by 
Schiissler & Vogler (2006) in magnetoconvection simulations 
of umbral dots. Interestingly, the parameter regime covered by 
those simulations is not the one relevant to sunspot penumbrae. 

Our aim here is to shed some light on the origin of dark-cored 
penumbral filaments and the surplus brightness of the penumbra. 
To that end we solve the 2D stationary heat transfer equation in 
a stratified uncombed penumbra formed by magnetic flux tubes 
in a stronger background field (Sects. [2] and [3]). The tubes carry 
an Evershed flow of hot plasma. Our calculations show that one 
such tube would be observed as a dark-cored filament due to 
the higher density of the plasma within the tube (Sect. |U). We 
also find that the Everhed flow heats the background atmosphere 
very efficiently, increasing its temperature. This suggests that 
the surplus brightness of the penumbra is due to the Evershed 
flow (Sect. 15.21) . Finally, we synthesize polarization maps us- 
ing the model atmospheres resulting from the simulations and 
compare them with polarimetric observations of dark-cored fila- 
ments (Sect. 15.31) . 



2. The model 

We describe a bright filament in the inner penumbra as a cylin- 
drical flux tube of radius R embedded in a stratified back- 



ground (umbral) atmosphere. The calculations are performed in 
a Cartesian coordinate system where the z-axis coincides with 
the vertical and the v-direction is defined such that the axis of 
the tube lies in the yz-plane. Since our primary goal is to identify 
the mechanism(s) responsible for the existence of dark cores, we 
adopt the simplest magnetic configuration possible, namely a po- 
tential field. The field is determined from the conditions V B = 
and V x B = in the xz-plane, i.e., we neglect variations along 
the tube axis because they are much smaller than variations per- 
pendicular to it (e.g., Borrero et al. 2004). 

In the tube's interior, the magnetic field is taken to be ho- 
mogeneous and directed along its axis, which is inclined by an 
angle y t to the vertical: B = B t (siny t e y + cosy^z)- Far from 
the tube we assume the background magnetic field to be of the 
form B = B^siny^ey + cosyb^X so that the field is homoge- 
neous with strength B\> and inclination y^. We use these condi- 
tions, together with the continuity of the radial components at 
the tube's boundary, to solve Laplace's equation in the xz-plane. 
Spectropolarimetry tells us that the field is weaker and more in- 
clined in the tubes (e.g., Bellot Rubio et al. 2004; Borrero et al. 
2004), hence we set B t < B^ and y t > TV The analytic expres- 
sions obtained from the calculations are given in Appendix A. 

Figure[T] shows the magnetic configuration of the model with 
R = 120 km, B t = 1200 G, y t = 90°, B h = 2000 G, andy b = 40°. 
The background field lines wrap around the tube (top left), lead- 
ing to enhanced field strengths on either side of the tube and 
reduced field strengths above and below it (top right). The in- 
clination of the background field also changes depending on the 
position. In particular, the field above the tube is always more 
horizontal than y\> due to the continuity of the normal compo- 
nents at the tube's boundary (bottom left). 
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We require lateral mechanical equilibrium, therefore the to- 
tal (gas plus magnetic) pressure on both sides of the interface 
separating the tube from the ambient medium is the same at the 
same geometrical height z. The magnetic field and the temper- 
ature then determine the gas density through the ideal gas law 
(with variable mean molecular weight to account for partial ion- 
ization). The equilibrium density distribution shown in the bot- 
tom right panel of Fig. [T] corresponds to the case in which the 
temperatures of the tube and the external medium are those dis- 
played in Fig. [2 As can be seen, a strong density enhancement 
occurs within the tube to compensate for its lower magnetic pres- 
sure. 

The condition of lateral force balance does not guarantee 
vertical equilibrium, so the tube may stretch in the vertical direc- 
tion. Using simple estimates, however, Borrero et al. (2006) have 
shown that the vertical stretching could be limited by buoyancy 
in the subadiabatic layers of the photosphere. Alternatively, non- 
potential fields may ensure vertical force equilibrium (Borrero 
2007). Since there is no generally accepted solution to this prob- 
lem, we restrict ourselves to the simple magnetic configuration 
described above in the hope that the results will not depend sig- 
nificantly on the exact topology of the field. 

3. Heat transfer equation 

To obtain the temperature distribution in the xz-plane we con- 
sider the stationary (d/dt - 0) heat transfer equation 



V-F = S, 



(1) 



neglecting gradients in the v-direction. Here, F = F r + F c rep- 
resents the heat flux vector and S the various energy sources, 
including Ohmic dissipation and the Evershed flow. The radia- 
tive flux F r is computed using the diffusion approximation 



f y = —k y vr 



(2) 



(e.g., Mihalas 1978), with k v the radiative thermal conductivity 
and T the temperature. Following Schlichenmaier et al. (1999), 
we take k v = 16 d~T 3 /(k^p), where & is Stefan-Boltzmann 
constant, &r the Rosseland mean opacity, p the gas density, 
and Dy the flux limiter originally introduced by Levermore & 
Pomraning (1981). The convective flux F c is evaluated using a 
linearized mixing length approach (Moreno-Insertis et al. 2002), 



F c = -k c [VT - (d77dz) ad ] = -Kc V c r, 



(3) 



where (dr/dz) a d = 7XdP/dz)V a d represents the adiabatic tem- 
perature gradient and V a( j the double-logarithmic isentropic tem- 
perature gradient. The computation of V a( j takes into account the 
local physical conditions and the partial ionization of hydrogen 
(Cox & Giuli 1968). In our simulations, V a( j varies between 0.1 
and 0.4 depending on the layer. The convective transport coeffi- 
cient k c is evaluated following Spruit (1977). To reduce the ef- 
ficiency of convection in the presence of magnetic fields we use 
r eff = 3888 K, z\ = 70 km, and k = 0.8 x 10" 11 . The calcula- 
tions are not very sensitive to the details of the energy transport 
by convection because k c is negligible in photospheric layers. 
We note, however, that a non-linear treatment of the convective 
flux, with k c varying with the changing superadiabaticity of the 
atmosphere, would be more appropriate to deal with large tem- 
perature perturbations (Sect. 14.21) . 

In the absence of penumbral flux tubes, the background at- 
mosphere is a plane-parallel stratified medium where S = and 
the physical parameters do not vary with x and y. Let us la- 
bel them with the subscript 0, i.e., 7q(z), k y q(z), k c o(z), and so 
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Fig. 2. Physical parameters in the cool umbral model of Collados 
et al. (1994). Clockwise from upper left, temperature, gas density, 
thermal coefficients, and energy fluxes. Radiative and convective 
parameters are shown in solid and dashed lines, respectively. 



on. Flux tubes embedded in this atmosphere present an obstacle 
to the heat flow because of their higher gas density, which de- 
creases k y . In addition, S ± within the tubes. As a result, the 
temperature gets modified from Tq(z) to T(x, z), but the gas pres- 
sure remains the same to ensure horizontal force balance. To find 
the new equilibrium configuration we solve Eq. (1) iteratively. 
Let Ti be the value of T(x, z) in the z-th step. This approximate 
solution does not satisfy Eq. (1) by an amount 



(4) 



where Si = S(Ti), K ti = K Y (Tt), and K ci = K c (Ti). Our goal is to 
find a perturbation 6T i+ \ such that the new temperature T i+ \ = 
Ti + STi+i leads to smaller errors, i.e., < e*. 

In terms of T i+ \, the heat transfer equation can be written as 



V • [-K Ti+1 VT i+1 - K ci+1 V c T i+ i] =S i+ i. 
Defining 

K TM = K r0 + 6oK Ti+ i, 
Kci+l = KcO + 8()Kci+U 



(5) 



(6) 



and neglecting the second-order terms £o%+iV(£7Vi) and 
SoK ci+ iV c (6Ti+i), Eq. ([5]) becomes 

v.[-^oV^ +1 -^ c0 v c ^ +1 ] = s i+ i-v-[-ic ri . + iVr I -iQ. I . +1 v c r f ].(7) 

The heat flux variation due to temperature variations is caused 
primarily by changes in the temperature gradient, while the ef- 
fects of changes in the transport coefficients are negligible to first 
order (Moreno-Insertis et al. 2002). This implies 



K l0 V6Ti+i » [K Ti+1 - K Yi ]VTi, 
K ci+l K ci 



(8) 
(9) 



allowing us to make the approximations /c n+1 ^ K Ti and K ci+ i ^ 
K ci m Eq. ©. With the additional assumption that Sj+i ^ Si, 
Eq. © can finally be rewritten as 



V • [-K r0 VdT i+ i - K c0 V c 6T i+ i] = 6i, 



(10) 



where the coefficients on the left-hand side depend only on z- We 
seek for periodic solutions of Eq. (ITOl) in the x-direction. Thus, 
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Fig. 3. Cut of the computational domain at x = km, showing 
a flux tube 120 km in radius. The vertical dotted lines mark the 
position of the 17 xz-planes where we solve the stationary heat 
transfer equation. The vertical solid line at y = 2083 km indi- 
cates the position of the xz-cut displayed in Figs. [Hand \5\ 



the following Neumann boundary conditions are imposed: 



5Ti(x, Z = Zmin) = Vx, 

6Ti(x = 0, z) = 6Ti(x = x max , z) = Vz, 



(11) 
(12) 



where the dots represent derivation with respect to x. The first 
condition prescribes zero temperature variations at the bottom 
of the computational domain, while the second forces the tem- 
perature perturbations to be minimum at the lateral boundaries. 
With these ingredients, Eq. dTOb is solved numerically for each 
height z using Fourier cosine transforms in the x-direction. 
Convergence (max|e;/Sj| < 10" 4 ) is achieved in 5 to 10 itera- 
tions. 



3.1. Simulation setup 

The computational domain is a cube extending 1000 km in the 
x-direction, 3100 km in the y-direction, and 1200 km in the z- 
direction. Initially, the cube is filled with an unperturbed umbral 
atmosphere (the cool model of Collados et al. 1994) with field 
strength = 2000 G and inclination = 40°. Figure 2 shows 
the run with depth of the temperature, gas density, and radiative 
and convective parameters in the unperturbed model. A magnetic 
flux tube is inserted in this atmosphere. It has a radius R = 120 
km, a current sheet 2 km thick, a field strength B t = 1200 G, 
and field inclinations varying between 45° and 87° as indicated 
in Fig. [3] The axis of the tube is placed at a height determined 
by the field inclination and the distance along the tube, starting 
with z = -326 km at y = km. 

The heat transfer equation is solved in 17 xz-planes at differ- 
ent y- values (Fig.0. Each plane is discretized in 501 x 601 grid 
points separated by 2 km. 



4. Results 

We first consider Ohmic dissipation as the only source of energy 
(Sect. l4~TT) . Thereafter, we examine how the results are modified 
by a hot Evershed flow along the magnetic flux tube (Sect. 14.21) . 



4.1. Ohmic dissipation as the only source of energy 

In this case, the source term of the heat transfer equation reduces 
to S = j 2 /cr, where j is the electric current associated with the 
jump of B at the tube's boundary (Eq. IA.4D and cr the electric 
conductivity, computed following Kopecky & Kuklin (1969). 

Figure 01 shows the simulation results at y = 2083 km. The 
tube blocks the energy coming from below and this produces a 
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Fig. 4. Top: temperature perturbations AT = T(x, z)-To(z) in the 
xz-plane at y = 2083 km. Ohmic dissipation is the only energy 
source considered. The circle represents the flux tube. Solid lines 
are lines of constant Rosseland optical depth. Bottom: contin- 
uum intensity at 487.8 nm emerging from the y = 2083 km 
plane, convolved with the Airy point- spread function of a 1-m 
telescope. The values are normalized to the quiet Sun continuum 
intensity Iqs- 



strong heating of its lower half and the adjacent medium. The 
associated temperature enhancements, however, are difficult to 
detect because they occur below the photosphere for the most 
part. 

The continuum intensity observed at the surface is deter- 
mined mainly by the temperature prevailing at Rosseland optical 
depth unity, tr = 1. In the top panel of Fig. |U lines of constant 
optical depth are indicated. As can be seen, the tr = 1 level is 
shifted upwards within the flux tube. The reason is the larger gas 
density in the tube compared with the external medium, which 
results from its weaker field strength and the condition of hori- 
zontal mechanical equilibrium. The enhanced density increases 
the opacity and moves the tr = 1 level to higher layers, where 
the temperature is lower owing to the stratification of the model. 
In accordance with the Eddington-Barbier relation, the contin- 
uum intensity decreases at the position of the tube. 

The bottom panel of Fig. |4] shows the continuum intensity 
emerging from the model at 487.8 nm, as computed by solv- 
ing the radiative transfer equation and convolving the result with 
the theoretical Airy point- spread function of a 1-m telescope. 
The flux tube is darker than its surroundings and do not ex- 
hibit bright edges, in clear contradiction with the observations. 
To reproduce the properties of dark-cored filaments, the walls of 
the tube should be hotter than the background atmosphere. The 
tr = 1 level would then encounter higher temperatures at the 
tube's boundary, generating two lateral brightenings. Ohmic dis- 
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Fig. 5. Same as Fig.ffl for the case in which the sources of energy 
are Ohmic dissipation and a hot Evershed flow. 

sipation in the current sheet, however, cannot provide the nec- 
essary heating because of the high conductivity of the plasma. 
In fact, one would have to reduce the conductivity by at least 
two orders of magnitude, or to consider extremely thin current 
sheets (a few meters thick), to heat the walls of the tube to the 
required degree. Even in that case the dark core would still show 
smaller continuum intensities than the surroundings. It is there- 
fore necessary to heat not only the walls of the tube, but also the 
tube itself. In this context, the Evershed flow represents a natural 
source of energy. 

4.2. Heating by the Evershed flow 

The Evershed flow is a radial outflow of mass associated with 
the more inclined magnetic fields of the penumbra (Title et al. 
1993; Stanchfield et al. 1997; Schlichenmaier & Schmidt 2000; 
Westendorp Plaza et al. 2001; Borrero et al. 2005; Bellot Rubio 
et al. 2006; Rimmele & Marino 2006; Ichimoto et al. 2007). 

An Evershed flow of hot plasma along a magnetic flux tube 
produces an energy flux Fe whose divergence can be evaluated 
from the entropy equation as 

V • F E = pcyVE [Vr - (d77<k)ad*zL (13) 

where cy is the specific heat at constant volume and Ve the flow 
velocity (ve = outside the tube). To account for the heating 
induced by the Evershed flow, we set S = j 2 /cr + V • F E in 
Eq. ©. We assume that the flow is parallel to the magnetic field 
(Bellot Rubio et al. 2004), and that its velocity changes with ra- 
dial distance as dictated by mass conservation due to the density 
decrease toward higher layers. 

The calculation of Eq. (IT3l) cannot be done in 2D because 
of the term V7\ To solve the problem we follow an iterative 
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Fig. 6. Section of the xz-plane at y = 2083 km showing gas den- 
sity perturbations Ap/p = [p(x,z) - Po(z)]/po(z) induced by 
the weaker field of the tube and the Evershed flow. Dash-dotted 
lines represent isotherms. Solid lines indicate lines of constant 
Rosseland optical depth. 

approach. V7 is estimated from the temperatures in the plane 
y = yj-i, T(x, yj-i,z), and a guess for the temperatures at y = yj, 
Tt(x, yj, z). Once S is known, the method described in Sect.[3]can 
be applied to the y = yj plane. This results in a new temperature 
distribution T i+ \(x, yj, z) which is used to update V7 for the next 
iteration. Convergence is reached in 10-50 steps. To initialize 
the calculations, the temperature distribution in the y = yj plane 
is taken to be the temperature at y = Vj_i, plus the temperature 
difference between the two planes in the non-Evershed case. 

The temperature excess induced by the flow at y = km is 
prescribed as a boundary condition. This condition determines 
the overall brightness of the flux tube and its surroundings, so it 
is relatively well constrained. In our calculations we use a tem- 
perature excess of 7500 K within the tube at y = km, to ensure 
that its brightness is compatible with the observations. An excess 
of 7500 K is also compatible with the maximum temperatures 
of -13 500 K that hot plasma rising adiabatically from the bot- 
tom of the convection zone would show at photospheric levels 
(Borrero 2007). 

The upper panel of Fig. \5\ displays the simulation results at 
y = 2083 km, assuming a flow velocity of ve = 7 km s" 1 in the 
first plane. The equilibrium temperature distribution is charac- 
terized by an intense heating of the tube and the surroundings, 
with temperature enhancements of up to 6000 K. This heating 
increases the brightness of the external medium and modifies 
the opacity of the plasma in a way that produces filaments with 
a dark core and two lateral brightenings (lower panel of Fig. [5]). 

Figure [6] examines in greater detail the origin of these struc- 
tures. We show the gas density perturbations, Ap/po = [p(x, z) - 
Po(z)]/po(z), in the xz-plane at y = 2083 km together with lines 
of constant temperature (dash-dotted) and constant Rosseland 
optical depth (solid). The heating caused by the Evershed flow 
shifts the tr lines to higher layers in and near the tube. For 
tr > 10, the lines of constant tr follow the isotherms closely, 
i.e., deep in the atmosphere the optical depth scale is primarily 
determined by the temperature, not by the gas density. This is 
due to the strong temperature dependence of the H" opacity. As 
the tube is approached from x = -150 km, the tr = 1 level 
crosses isotherms corresponding to higher temperatures, which 
explains the larger continuum intensities of the filament com- 
pared with its surroundings. Very close to the tube's boundary 
the tr = 1 line encounters regions of even higher tempera- 
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Fig. 7. Dark-cored penumbral filament produced by the flux 
tube considered in this work as seen in continuum intensity at 
487.8 nm. The Evershed flow has a velocity of ve = 7 km s" 1 
and induces a temperature excess of 7500 K within the tube at 
y = km. 



tures, producing a maximum in the continuum intensity (the lat- 
eral brightenings). As soon as the tr = 1 line enters the tube it 
evolves in a region of enhanced gas density. The increased opac- 
ity moves the tr = 1 surface upward by about 40 km, where the 
temperature is cooler. This generates a dark core. 

The density enhancement is more prominent near the top of 
the tube because the difference of magnetic pressures between 
the interior and the external medium attains a maximum there 
with respect to the external gas pressure. The temperature per- 
turbations induced by the Evershed flow decrease the large den- 
sity contrasts resulting from the weaker field of the tube (bottom 
right panel of Fig. Q]), but this effect is less important in deter- 
mining the equilibrium density distribution. 

Two interesting facts deserve further consideration. First, the 
tr = 0.1 level is reached at the top of the tube, some 40 km 
above tr = 1. Thus, a substantial fraction of the line forma- 
tion region is contained within the flow channel. Second, the 
density distribution depicted in Fig. [6| tends to stabilize the tube 
against vertical stretching, since it causes negative buoyancy at 
the top of the tube. A similar density distribution has been found 
by Borrero (2007) under the assumption of flux tubes in exact 
force balance. 

Figure [7] displays a continuum (487.8 nm) image of the 
tube as it would appear through a 1-m telescope at disk center. 
Clearly, the morphological properties of dark-cored filaments 
are well reproduced: the inner footpoint shows up as a bright 
penumbral grain, the central dark lane is surrounded by two lat- 
eral brightenings separated by roughly 250 km (~073), the in- 
tensity ratio between the dark core and the lateral brightenings 
is 0.93, and the dark core can be traced for more than 2000 km 
(Scharmer et al. 2002; Sutterlin et al. 2004; Rouppe van der 
Voort et al. 2004; Bellot Rubio et al. 2005; Rimmele 2008). We 
also note that the distance between the lateral brightenings re- 
flects the true diameter of the underlying magnetic flux tube, at 
least at the resolution of a 1-m telescope. 



5. Discussion 

5.1. Origin of dark-cored penumbral filaments 

Schlichenmaier et al. (1999) studied the radiative cooling of 
hot flux tubes surrounded by an initially isothermal atmosphere. 
They found that the tubes cool down quickly in the absence of 
energy sources, reaching thermal equilibrium with the external 
medium in only a few tens of seconds. In the case of optically 
thick tubes (initial temperatures above 10000 K), a steep cool- 
ing front develops and migrates towards the tube axis at constant 
velocity, while optically thin tubes (initial temperatures below 
7500 K) cool more or less homogeneously over their cross sec- 
tions. In either case, no dark cores would be observed. 



Our simulations can be regarded as an extension of the 
work by Schlichenmaier et al. (1999). The main difference is 
that we consider a stratified atmosphere consisting of magnetic 
flux tubes in a stronger ambient field. Under these conditions, 
the equilibrium temperature distribution is no longer symmetric 
around the tube axis. We have demonstrated that one such tube 
would exhibit a central dark core due to the higher density (hence 
larger opacity) of the plasma within the tube, which moves the 
tr = 1 level to higher (cooler) layers. The same mechanism pro- 
duces dark lanes in umbral dots (Schtissler & Vogler 2006) and 
gappy penumbrae (Spruit & Scharmer 2006), although in these 
structures the density enhancement is a consequence of overturn- 
ing convection. In an uncombed penumbra, the density increase 
is essentially due to the weaker field of the tubes. 

It is important to remark, however, that the observations can- 
not be fully explained without a hot Evershed flow along the 
tubes. The flow increases the brightness of the flux tube relative 
to the surroundings and generates the two lateral brightenings 
of the filaments; in the absence of this energy source, the tubes 
would actually be darker than their environs. 

Borrero (2007) has also been able to reproduce dark-cored 
filaments using thick penumbral flux tubes in magnetohydro- 
static equilibrium. The condition of vertical and horizontal force 
balance determines the temperature distribution in and around 
the tubes, given the magnetic field distribution. This results in 
structures that are hotter than the external medium at the same 
geometrical height, in good qualitative agreement with our cal- 
culations. Borrero (2007) did not identify the source of the heat- 
ing, but we hypothesize that the higher temperatures required 
to maintain the tubes in force equilibrium are actually produced 
by the Evershed flow. To verify this conjecture, the energy and 
momentum equations must be solved simultaneously in the pres- 
ence of hot upflows. 

5.2. Surplus brightness of the penumbra 

An important result of the calculations described in Sect. 14.21 is 
that the background atmosphere itself is much brighter than it 
would be in the absence of an Evershed flow. The continuum 
intensity without flows is only 0.06 /qs (Fig. [4]), corresponding 
to a very cool umbra. In contrast, when a hot upflow is consid- 
ered, the background shows intensities of up to -0.7 Iqs near the 
tube (Fig. [7]). This enormous difference can explain the surplus 
brightness of the penumbra relative to the umbra. 

The average continuum intensity emerging from the box of 
Fig.[7]is ~0.5 /qs, i.e., slightly smaller than the observed penum- 
bral brightness but of the same order of magnitude. Adjustments 
of the temperature in the background atmosphere, the geom- 
etry of the flux tubes, and/or the boundary conditions for the 
Evershed flow should easily lead to a closer match. 

5.3. Polarimetric signatures of dark-cored penumbral 
filaments 

In this section we investigate the appearance of penumbral flux 
tubes in polarized light using the results of Sect. 14.21 We have 
chosen the Fei line at 630.25 nm to solve the radiative transfer 
equation because many instruments, including the SOUP mag- 
netograph at the Swedish Solar Telescope (SST) and the po- 
larimeters of the Solar Optical Telescope aboard Hinode, have 
measured the properties of dark-cored filaments in this line. 

The four top panels of Fig. [8] show the flux tube as it would 
be recorded through a 1-m telescope at disk center in continuum 
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Fig. 8. Top: polarization maps of the dark-cored penumbral fil- 
ament resulting from the simulation of Sect. 4.2 in the Fei 
630.25 nm line. The four panels show the continuum intensity 
and the wavelength-integrated (unsigned) Q, U, and V polariza- 
tion through a 1-m telescope at disk center. The polarization pro- 
files have been normalized to the quiet Sun continuum. Bottom: 
magnetogram at +10 pm from line center and Dopplergram at 
±15 pm from line center. 



intensity and wavelength-integrated Stokes Q, U, and V polar- 
ization, respectively. These maps can be compared with real ob- 
servations. The dark lane is more prominent in polarized light 
(2, U, and V) than in continuum intensity, as observed with 
Hinode (Bellot Rubio et al. 2007) and the SST (van Noort & 
Rouppe van der Voort 2008). In total circular polarization, how- 
ever, the lateral brightenings disappear at too short a distance 
from the filament head. We believe this is due to the geome- 
try assumed for the flux tube, which never becomes horizontal. 
Because of the heights attained by the tube, the radiative cooling 
of the Evershed flow proceeds at a very fast pace. Also, the tube 
leaves the region of maximum sensitivity of Stokes V compara- 
tively soon, which decreases the circular polarization signal. For 
these reasons, we expect that a better choice of the tube incli- 
nation at large radial distances will produce longer filaments in 
circular polarization. 

The bottom panels of Fig.[8]display a synthetic magnetogram 
and a Dopplergram of the filament. They have been constructed 
using the Stokes V profile of Fei 630.25 nm at +10 pm from 



line center, and the corresponding Stokes / profile at +15 pm, 
respectively. The magnetogram signal is weaker in the dark core 
than in the lateral brightenings (white means less signal), sug- 
gesting weaker and/or more inclined fields. This is exactly what 
has been inferred from polarimetric measurements of dark-cored 
filaments (Langhans et al. 2007; Bellot Rubio et al. 2007; van 
Noort & Rouppe van der Voort 2008). The Dopplergram shows 
a strong signal in the dark core, indicating the existence of a 
flow there. Given the disk-center position of the tube and the 
nearly horizontal inclination of the flow, the Doppler shift is not 
very conspicuous in this particular example; observations closer 
to the limb would certainly show stronger signals. In any case, 
the Dopplergram of Fig. [8] provides the same information as 
the high-resolution spectroscopic measurements of Bellot Rubio 
et al. (2005), Rimmele & Marino (2006), and Langhans et al. 
(2007): the flow appears preferentially in the dark core, not in 
the lateral brightenings or the adjacent medium. 

The qualitative agreement between synthetic and observed 
parameters is not surprising. As mentioned in Sect. 14.21 a signif- 
icant fraction of the line-forming region is contained within the 
flux tube. This means that the physical conditions of the tubes 
leave their signatures in the Stokes spectra of magnetically sensi- 
tive lines. Since the tubes have weaker and more inclined fields, 
together with strong flows, the polarization profiles emerging 
from them cannot indicate otherwise. To some extent these ar- 
guments justify the conclusions drawn from simple interpreta- 
tions of polarimetric measurements, but we caution that precise 
determinations of the magnetic field vector and flow velocity of 
dark-cored penumbral filaments will require sophisticated inver- 
sion techniques such as those used by Borrero et al. (2005) or 
Jurcak et al. (2007). 

6. Conclusions 

The heat transfer and radiative transfer calculations presented in 
this paper support the concept of a penumbra formed by small 
(but optically thick) magnetic flux tubes that carry hot flows, 
as deduced from high-resolution observations and spectropolari- 
metric measurements (see Solanki 2003 and Bellot Rubio 2004 
for reviews). Tubes about 250 km in diameter explain not only 
the existence of dark-cored filaments, but also the surplus bright- 
ness of the penumbra; the Evershed flow efficiently heats the 
plasma outside the tubes, increasing its temperature to values 
compatible with the observations. 

Further improvements of the model should include a more 
realistic treatment of the magnetic topology of the tubes and the 
external atmosphere, a better description of convection, and a 
full 3D solution of the heat transfer equation. In our opinion, 
however, these improvements will not change the conclusion that 
the uncombed model is the best representation of the penumbra 
at our disposal. 
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Appendix A: Magnetic field distribution 

We take advantage of the information provided by observations 
and numerical simulations to constrain the magnetic properties 
of the model. As mentioned in Sect. 2, the field strength in- 
side the tube, B u must be smaller than in the external medium, 
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#b- Also, the inclination of the field in the tube, y t , should be 
larger than in the background far from it, yt,. These inclinations 
are measured with respect to the local vertical, as sketched in 
Fig. [AH 

A potential field satisfies the conditions V • B = and V x 
B = 0. We solve this problem in cylindrical coordinates rOy' , 
where the /-axis coincides with the axis of the tube (Fig. IA.lD . 
In the tube's interior (r < R - S), the field is assumed to be 
homogeneous and directed along its axis, i.e., 



By 



Be 
B t . 



0, 



(A.l) 



Variations of the field along the / -axis are ignored because they 
are much smaller than variations perpendicular to the tube's axis 
(Borrero et al. 2004). In the background atmosphere, we assume 
that far from the tube the magnetic field vector is inclined by an 
angle y' h = 90 + yb - y t with respect to the z'-axis, and that the 
field does not depend on z' . Further, we require continuity of the 
radial component of the field at r = R. Using these boundary 
conditions, the potential solution for the magnetic field in the 
external medium (r > R) reads 



B r = B b cosy b sin 0(1 -R 2 /r 2 ), 
B e = B h cosy; cos<9(l + R 2 /r 2 ), 
By = B h siny b 

Finally, in the current sheet (r e [R - 
By to vary linearly, so that 



(A.2) 



S,R]) we demand B e and 



B r = 0, 

B e = 2B h cosy b cos# [r - (R - 6)] /6, 
By = B t + (B b siny b - B t ) [(r - (R - S)] 1 6. 



(A.3) 



With this choice, the magnetic field vector is continuous across 
the walls of the tube and verifies the condition V • B = 0. 

The currents associated with this magnetic field are non-zero 
only in r e [R - S, R] and can be expressed as 



h 



0, 



Je = — 



B h siny b - B t 



(A.4) 



Jy 



2B b 



R-S 



cos y b cos 0, 



li being the magnetic permeability in vacuum. The current in 
the direction of the tube axis, jy , is proportional to the jump of 
the vertical component of the magnetic field vector. It reaches a 
maximum at z' =0 and vanishes at the top and bottom of the 
tube (where = 90° and = 270°). 
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